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Abstract 


At  the  request  of  the  Illuminating  Engineering  Research  Institute 
of  the  Illuminating  Engineering  Society,  the  National  Bureau  of  Standards 
has  written  a computer  program,  for  use  on  the  IBM  650,  high  speed  digital, 
computer,  which  converts  into  terms  of  the  1931  CIE  (International 
Commission  on  Illumination)  standard  observer  and  coordinate  system  data  of 
spectral  transmittance  or  spectral  directional  reflectance  of  a sanple 
illuminated  by  any  light  source,  including  those  sources,  such  as  fluores- 
cent lamps,  which  have  spectral  lines  superimposed  on  the  continuous  spectrum 
of  the  sourceo  A description  of  the  data  required  by  the  program  is  given 
together  with  detailed  instructions  for  coding  the  data. 


- 2 - 


i 


Conten  ts 


I«  General  description  of  the  computing  program 
II*  Form  of  the  input  data 

III.  Coding  the  data 

IV.  Operation  of  the  IBM  6^0  computer 

V.  Form  of  the  output  data 

VI.  Summary 

VII .  Ackn  owledgeme  nts 

VIII.  B ibli 0 gr aphy 
Appendix  A 
Appendix  B 
Appendix  C 
Appendix  D 
Appendix  E 


Page 

3 

h 

6 

9 

10 

10 

10 

10 

12 

20 

22 

25 


I 


. ^ ^ 


lo  General  description  of  the  conputing  program. 

The  purpose  of  this  computing  program  is  the  conversion  of  spectral 
directional  reflectance  or  spectral  transmittance  data  into  colorimetric 
terms  of  the  1931  CIE  (international  Commission  on  Illumination)  standard 
observer  and  coordinate  system  The  wavelength  range  covered  is 

380  to  770  millimicrons.  Any  light  source  may  be  used  in  the  confutations 
including  those  sources^  such  as  fluorescent  lamps,  which  have  spectral 
lines  superimposed  on  the  continuous  spectrum  of  the  source. 

Up  to  $ spectral  lines  can  be  included  in  the  data  for  the  source.  It 
is  necessary  to  give  only  the  wavelengths  at  which  the  spectral  lines  occur 
and  the  relative  radiant  fluxes  of  these  lines.  The  computer  program  per- 
forms a third-difference,  osculatory  interpolation  on  the  data  of  the 
CIE  standard  observer  and  the  spectral  directional  reflectance  or  spectral 
transmittance  data  of  the  sample  at  the  proper  wavelengths.  The  third- 
difference,  osculatory  interpolation  formula  has  the  condition,  however, 
that  no  such  interpolation  can  be  made  within  the  first  and  last  wavelength 
interval,  namely  38O  to  390  millimicrons  and  760  to  770  millimicrons.  Since 
the  possibility  exists  that  a spectral  line  could  fall  within  the  first  or 
last  wavelength  interval,  the  program  has  been  written  to  perform  linear 
interpolation  in  these  wavelength  regions. 

T«ro  possibilities  for  handling  data  are  available.  The  conversion  to 
CIE  colorimetric  terms  can  be  made  for  (1)  a series  of  source-sample  combina 
tions  where  neither  the  data  of  the  source  nor  the  data  of  the  sample  are 
repeated,  or  (2)  a number  of  sanfles  relative  to  the  same  source. 

The  resulting  data  of  the  computations  are  punched  on  a single  card 
giving  the  identification,  tristimiiLus  values,  and  chroraaticity  coordinates 
of  the  source-sample  combination. 

A complete  listing  of  the  FOR  TRANSIT  [3l  code,  used  to  assemble  the 
program  deck,  is  shown  in  Appendix  A.  Since  the  FOR  TRANSIT  system,  used 
on  the  IBM  6^0  computer,  is  coirpatible  with  the  FORTRAN  j^Ii, ^3  system,  used 
on  the  IBM  7OI1  computer,  it  was  possible  to  prepare  the  computing  program 
on  the  IBM  70I4  computer  at  the  National  Bureau  of  Standards,  and  with  minor 
modifications,  obtain  a program  which  could  be  assembled  on  the  IBM  6^0 
computer  at  the  U.  S.  Department  of  Agriculture.  Thus,  the  FORTRAN  program, 
described  in  NBS  Report  6613,  performs  the  same  computations  as  the  FOR 
TRANSIT  program  described  herein.  The  modified  FORTRAN  [63  program  which 
was  assembled  on  the  IBM  70k  computer  is  shown  in  Appendix  B. 

The  program  deck,  supplied  with  these  instructions,  is  this  FOR  TEIANSIT 
code  assembled  in  "machine  language"  on  a basic  IBM  650  computer  equipped 


^ Numbers  in  brackets  refer  to  bibliography  on  page  10  of  this  report. 
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with  a special  character  device  (Group  II),  using  the  FOR  TRANSIT  l(S)  [3^ 
compiling  program,  and  the  FOR  TRANSIT  ^33  control  panel. 

The  computing  program  has  been  checked  against  data  computed  manually 
on  desk  calculators.  The  results  of  the  two  sets  of  calculations  agree 
exactly  to  at  least  five  decimal  places  which  is  considered  sufficient  for 
the  problems  to  be  computed  by  this  program. 

The  running  time  on  the  computer,  after  the  program  deck  has  been 
loaded,  is  approximately  1 minute  1^  seconds  per  source-sample  combination. 

II.  Form  of  the  input  data. 

The  input  data  are  punched  on  standard  80  column  IBM  cards  and  arranged  1 
in  the  following  order. 

(1)  Control  card, 

(2)  Continuum  of  the  source. 

(3)  Number  of  spectral  lines  to  be  considered. 

(h)  Wavelengths  of  the  spectral  lines. 

(5)  Radiant  fluxes  of  the  spectral  lines. 

(6)  Spectral  directional  reflectance  or 
spectral  transmittance  of  sample. 

The  information  listed  above  are  punched  in  the  first  70  columns  of  the  ^ 
cards  in  accord  with  instructions  to  be  described  in  detail  later  in  this 
section.  Columns  71  through  80  are  not  read  by  the  computer  and  may  be  used 
to  number  the  cards  sequentially,  or  for  other  identifying  information. 

In  Appendix  D is  shown  the  listing  of  the  deck  of  input  data  cards  used 
to  check  this  program.  Note  that  the  data  shown  on  any  line  in  Appendix  D 
are  punched  on  a single  card.  The  column  numbers  are  shown  for  reference 
only  as  the  first  two  lines  at  the  top  of  each  page.  These  column  number 
cards  do  not  appear  in  the  original  input  data  deck  used  for  the  computations. 
Each  card  has  been  given  a number  in  columns  7U  and  75 • These  card  numbers 
will  be  used  in  the  description  which  follows  for  each  type  of  input  data 
format. 

The  IBM  650  uses  a fixed  input  and  output  format  consisting  of  eight 
10  column  fields  per  card.  A field  of  a punch  card  is  a specific  number  of 
columns  in  which  data  may  be  placed.  Thus,  the  first  field  consists  of 
columns  1 through  10,  the  second  field  columns  11  through  20,  and  so  forth 
until  the  eighth  field  is  reached  which  consists  of  columns  71  through  80. 

(1)  Control  card. 

The  control  card,  for  example  cards  numbers  1,  17,  33,  U9,  56,  and  63, 
Appendix  D,  serves  to  control  transfers  within  the  program,  depending  upon 
the  input  data  to  be  used,  and  also  to  identify  the  input  and  output  data 


with  a particular  source  “Sample  combination. 

The  first  field  of  the  control  card,  columns  1 through  10,  is  punched 
with  "0"  in  columrsl  through  9 and  with  a "1",  "2”,  or  ”3“  in  column  10.  The 
punch  in  column  10  serves  to  control  the  transfers  within  the  program  and 
has  the  following  significance. 


Punch  in 

Column  10  Significance 


1 Read  in  data  of  new  source  and 
data  of  new  sau^le. 

2 Read  in  data  of  new  sample  but  use 
the  data  of  the  last  given  source. 

3 Stop. 


The  source-sample  combination  is  identified  by  a numerical  designation 
punched  in  the  second  field,  the  last  digit  appearing  in  column  20  and  all 
unused  columns  being  punched  ”0".  The  second  field,  on  output,  is  punched 
on  each  answer  card  and  is  therefore  needed  to  correlate  the  output  data 
with  the  source“sample  combination. 

(2)  Data  for  the  continuum  of  the  source. 


The  spectral“irradiance  data  of  the  continuum  of  the  light  source,  at 
10  millimicron  intervals  over  the  wavelength  range  380  to  770  millimicrons, 
are  punched  in  the  first  seven  fields  of  the  card  with  the  exception  of  the 
last  card  which  has  the  first  five  fields  punched.  See  cards  nuiribeis  2 
through  7>  Appendix  D.  There  must  be  i^O  values  of  spectral  irradiance  of 
the  continuum  given,  one  value  for  each  of  the  1|0  wavelengths  from  38O  to 
770  millimicrons,  inclusive.  If  data  of  spectral  irradiance  are  not  avail- 
able  for  any  particular  wavelengths,  it  is  necessary  to  extrapolate  data 
for  these  wavelengths.  If  the  extrapolated  data  are  very  small  in  magnitude, 
zero  value  may  be  assigned  to  the  spectral  irradiance  at  those  wavelengths. 
This  can  be  done  by  punching  "0"  in  all  ten  columns  of  that  field.  There 
must  be  six  cards  for  the  spectral  “irradiance  data  on  the  continuum. 

(3)  Kuniber  of  spectral  lines  to  be  used  in  the  computation. 

This  card,  for  example  card  number  8,  Appendix  D,  tells  the  computer 
the  number  of  spectral  lines  to  be  used  in  the  computation.  Up  to  $ spectral 
lines  may  be  used.  At  least  one  spectral  line  must  be  used.  The  card  is 
punched  in  the  first  field,  columns  1 through  10.  A ‘’0”  is  punched  in 
columns  1 through  9 and  the  number  of  lines  is  punched  in  colium  10. 

(U)  Wavelengths  of  the  spectral  lines  present  in  the  source. 


The  wavelengths,  in  millimicrons,  of  the  spectral  lines  present  in  the 
source  are  punched,  as  in  card  number  9>  Appendix  D,  in  the  first  five  fields 


of  a single  card.  If  less  than  five  wavelengths  are  present,  it  is  necessary 
to  punch  "0"  in  all  ten  columns  of  the  unused  fields. 


{$)  Radiant  flux  data  of  the  spectral  lines. 


After  the  wavelengths  of  all  of  the  spectral  lines  present  in  the  source 
have  been  read  into  the  computer,  -the  radiant  flux  of  each  line  relative  to 
the  continuum  is  read  in.  The  relative  radiant  flux  of  the  spectral  line  is 
defined  as  the  peak  height  of  the  spec troradiome trie  curve  for  the  spectral 
line,  above  the  continuum  of  the  source,  deterndned  for  a slit-width  of  10 
millimicrons,  the  wavelength  interval  used  for  summation  in  this  program. 
These  data  are  punched  in  a single  card  as  in  card  number  10,  Appendix  D,  in 
the  first  five  fields  of  the  card.  If  less  than  five  wavelengths  are  present, 
it  is  necessary  to  punch  "0"  in  all  ten  columns  of  the  unused  fields.  There 
must  be  the  same  number  of  data  of  relative  radiant  flux  as  there  are  wave- 
lengths given,  and  the  flux  data  must  be  punched  in  the  same  order  as  the 
wavelength  data. 


If  it  is  desired  to  make  the  colorimetric  conversion,  ignoring  the  con- 
tribution of  the  spectral  lines,  it  is  necessary  only  to  assign  zero-flux 
values  to  all  of  the  lines.  The  condition  stated  in  (3)  above,  that  at 
least  one  line  must  be  given,  still  holds.  Thus,  a single  line  must  be  in- 
diqated  with  its  wavelength,  but  with  a zero  flux. 


(6)  Spectral  data  of  the  sample. 


The  spectral  directional  reflectance  or  spectral  transmittance  data  of 
the  sample,  at  10-millimicron  intervals  over  the  wavelength  range  380  to  770 
millimicrons,  are  punched  in  the  first  seven  fields  of  the  card  with  the  ex- 
ception of  the  last  card  which  has  the  first  five  fields  punched.  Se^  for 
exanple,  cardsnuiribers  11  through  16,  Appendix  D.  There  must  be  hO  values  of 
spectral  directional  reflectance  or  spectral  transmittance  given.  If  data 
are  not  available  for  a^y  wavelengths,  for  example  380,  390,  760,  or  770 
millimicrons,  it  is  necessary  to  extrapolate  values  for  these  wavelengths. 
There  must  be  six  cards  of  spectral  data  for  the  sample. 


III.  Coding  the  data. 


The  above  discussion  indicates  the  way  in  which  the  data  are  to  be  pre- 
pared on  punched  cards  for  introduction  into  the  computer.  This  section  will 
deal  with  the  coding  of  -the  data  by  the  scientist  or  engineer  for  use  by  the 
key-punch  operator  in  preparing  the  cards.  Since  the  key-punch  operator  is 
usually  located  at  the  computing  center,  away  from  the  laboratory,  and  has 
no  previous  knowledge  of  the  problem,  it  is  necessary  to  write  out  the  data 
of  the  problem  in  a form  that  the  key-punch  operator  can  follow  with  the 
least  chance  of  making  errors  in  punching. 


It  is  suggested  that  the  data  be  written  on  squared  paper  or  graph  paper 
which  is  at  least  80  squares  wide.  The  squares  at  the  top  of  the  sheet 
should  be  numbered  from  1 through  80  and  the  eight  fields,  10  columns  wide,|pi 


“ 7 “ 


should  be  indicated  by  -vertical  rulings » A ”+"  (plus  or  *'12"  punch)  must 
be  punched  in  column  73»  This  can  be  accomplished  by  gang-punching  the 
in  all  cards  to  be  used  in  the  input  data  deck  before  they  are  gi-ven  to  the 
key-punch  operator,  or  by  ha-ving  the  key-punch  operator  overpunch  (put  a 
”12’'  punch  in)  column  73  when  the  whole  card  is  punchedo 


The  arrangement  of  the  data  on  cards  has  been  covered  in  II  above  for 
the  control  card  and  the  card  indicating  the  number  of  spectral  lines.  For 

the  remaining  cards,  the  data  must  be  entered  in  the  following  manner.  Data 

are  located  in  the  first  seven  fields  of  the  card.  If  negative  -values  are 

needed  for  the  data,  the  minus  sign  is  indicated  by  an  "11"  overpunch  in  the 

units  column  of  the  field  containing  the  negative  value,  for  example,  col- 
uims  10,  20,  30,  and  so  forth.  All  data  must  be  represented  as  decimal  num- 
bers in  the  form  .xxxxxxxxPP,  where  "x"  represents  the  da-ta  and  the  "PP"  is 
$0  plus  the  power  of  10  which  is  needed  to  shift  the  decimal  point  back  to 
its  proper  position  in  the  original  data.  The  decimal  poin-fcs  are  never 
punched.  All  blank  columns  must  be  punched  with  "0"  between  columns  1 through 
70o  Thus,  the  following  examples  would  be  coded  as  shown. 


Original 

Data 

loli6 

,lii6 

,01i|6 


Equi -valent 
Form 

olh6  X 10^ 
olli6  X lOO 
„lii6  X 10“1 


Coded 

Data 

lii6000005l 

lii600000^0 

lU600000li9 


For  coding  purposes,  the  original  data  on  the  continuum  of  the  source 
and  the  fluxes  of  the  spectral  lines  are  to  be  considered  as  numbers  on  a 
scale  whose  maximum  is  about  100,  that  is,  96,3  and  not  0.963;  the  original 
wa-velength  data  are  to  be  considered  as  always  being  in  millimicrons,  that 
is,  for  example  ii0U,7  (or  any  number  of  decimal  places  up  to  five);  and  the 
original  spectral  data  of  the  sample  are  to  be  considered  as  always  being 
in  terms  of  decimal  numbers,  such  as  0,17h2, 

In  addition  to  the  data  of  the  source-san^le  combinations  to  be  computed, 
it  is  necessary  to  read  into  the  computer  the  data  of  the  CIE  standard  ob- 
server and  the  wavelength  range,  380  to  770  millimicrons,  over  which  the 
computations  are  to  be  made.  In  Appendix  C are  shown  the  2k  cards  which  com- 
prise this  da ta-cons tacts  deck,  Car<^ numbers  1 through  6 contain  the  da-ba 
of  the  X tristimulus  function;  7 through  12,  the  data  of  the  Y -tristimulus 
function;  13  through  18,  the  data  of  the  Z tristimulus  function;  and  19 
through  2J4,  the  wavelengths  for  each  10-millimicron  inter-val  between  38O 
and  770  millimicrons.  These  2h  cards  may  be  made  a part  of  the  program  deck 
or  be  considered  as  a part  of  the  input  data.  In  ei-ther  case,  care  should 
be  taken  to  see  that  the  data  constants  are  included  only  once  for  any  pass 
through  the  coirputer. 

The  arrangement  or  order  of  the  cards  must  now  be  considered  since  the 
cards  must  be  in  the  proper  order  -to  insure  operation  of  the  program,  A 
generalized  discussion  follows  of  the  two  types  of  data  handling  possibilities 
available  with  this  program  and  how  the  cards  are  arranged. 


The  first  possibility  available  is  to  read  into  the  computer  a new 
source  and  a new  sample  for  each  computation.  ¥e  shall  call  this  Case  1. 

The  computer  is  informed  that  Case  1 type  of  data  input  is  to  be  used  by  a 
*'l"  punch  in  column  10  of  the  control  card.  The  data  deck  for  Case  1 would, 
in  general,  look  like  the  following : 

Control  Card  OOOOOOOOOl  (in  first, field) 


Continuum  data  (6  cards) 

Number  of  spectral  lines  in  source  (1  card) 

Wavelengths  of  the  spectral  lines  (1  card) 

Fluxes  of  the  spectral  lines  (1  card) 

Spectral  data  of  the  sample  (6  cards). 


The  data  deck  of  the  type  for  Case  1 is  always  used  for  the  first  source- 
sample  combination  of  a particular  series  of  computations  since  it  is  the 
only  possibility  which  allows  the  computer  to  read  in  both  a new  source  and 
a new  sample. 

The  second  possibility  available  is  to  read  into  the  computer  only  the 
data  for  a new  sample  and  perform  the  computation  using  the  source  which  was 
last  read  in.  We  shall  call  this  Case  2.  A ”2"  punch  in  column  10  tells  the 
computer  to  save  in  memory  the  last  source  read  into  the  computer  and  use  it 
with  the  new  sanple  data  which  will  be  read  in  to  perform  this  computation. 
The  data  deck  for  Case  2 would,  in  general,  look  like  the  following; 

Control  Card  0000000002  (in  firs  t field) 

Spectral  data  on  sample  (6  cards). 


After  all  of  the  data  to  be  computed  have  been  read  into  the  computer 
and  the  computations  performed,  it  is  necessary  to  inform  the  computer  that 
the  problem  is  finished.  This  is  done  by  placing  a card,  with  “0000000003" 
punched  in  the  first  field,  and  " 0000000000"  in  the  second  field,  after  the 
last  data  card  in  the  data  deck.  This  we  shall  call  Case  3«  When  this  card 
is  reached  and  read  into  the  computer,  a transfer  is  made  to  a stop  instruc- 
tion in  the  program  and  no  further  instructions  are  available  to  the  computer. 


Thus,  a complete  deck  might,  in  general  terms,  look  something  like  the 
following: 


Program  Deck 

X (read  in  X tris timulus  function) 

I (read  in  Y tri stimulus  f unc  tion ) 

Z (read  in  Z tristimulus  function) 

Lambda  (read  in  wavelengths  380,  390,  ...,  ?60,  770) 

Case  1 (read  in  source  M and  sample  A) 

Case  1 (read  in  source  N and  sample  B) 

Case  2 (read  in  sample  C,  use  source  N) 

Case  2 (read  in  sample  D,  use  source  N) 

Case  2 (read  in  sample  E,  use  source  N) 

Case  1 (read  in  source  0 and  sample  F) 

Case  2 (read  in  sample  G,  use  source  0) 

Case  3 (stop). 
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The  data  used  to  check  this  programj  listed  in  Appendix  D,  show  how  an 
actual  problem  was  set  up  for  running  in  the  computer.  The  problem  in 
Appendix  D is  really  $ sub-problems  which  are  identified  in  columns  19  and 
20  of  the  control  cards  as  sub-problems  9 through  13-  Sub-problem  9 covers 
cards  numbers  1 through  l6.  This  is  a Case  1 type  of  data  input  (note  the 
"1”  punch  in  column  10  of  card  number  1)  introducing  a new  source  and  a new 
sample.  Sub-problems  10  and  11  cover  cards  numbers  1?  through  32  and  33 
through  U85  respectively.  The  data  input  for  both  of  these  is  of  the  Case  1 
type.  Sub-problems  12  and  13,  cards  numbers  U9  through  $$  and  56  through  62, 
respectively,  are  recomputations  of  sub-problems  10  and  11  respectively  but 
of  the  Case  2 type  of  data  input  using  the  data  of  the  source  read  in  for 
sub-problem  11.  Sub-problem  13  is  the  last  conversion  to  be  made  during  this 
run  on  the  computer^  and  therefore,  following  the  last  data  card  of  sub- 
problem 13,  which  is  card  number  62,  we  have  the  Case  3 card  (card  number  63) 
which  tells  the  computer  that  it  has  reached  the  last  card  of  this  computing 
problem  and  that  it  is  to  stop. 

IV.  Operation  of  IBM  650  computer. 

The  following  sequence  of  operations  will  cause  the  coitputer  to  execute 
the  complete  conputation,  provided  no  errors  are  present  in  the  data  deck. 

1)  Set  the  console  as  follows? 


2)  Insert  FOR  TRANSIT  533  control  panel. 

3)  Load  blank  cards  in  punch  hopper. 

i;)  Set  Storage  Entry  Sign  Switch  to  PLUS. 

5)  Load  program  deck  into  read  hopper. 

6)  Load  input  data  deck  into  read  hopper  arranged  in  the  order  given 
in  III.  Coding  the  data,  above,  behind  the  program  deck  and  add 
three  blank  cards  afier  the  deck. 

7)  Press  "RESET"  key  on  console. 

8)  Press  "PROGRAM  START"  key  on  console. 

The  computer  reads  in  the  program  deck,  the  CIE  tristim\iLus  function  data, 
the  wavelengths,  and  the  first  source-sample  combination,  performs  the  com- 
putation, punches  the  answer  card,  then  reads  the  second  source-sample  com- 
bination, and  so  forth  until  the  last  card  is  reached. 


Storage  Entry?  70  1952  9999 
Switches?  Programmed 


STOP 

RUN 

RUN 

UPPER 

SENSE 

STOP 


Half  Cycle 


Control 

Display 

Overflow 

Error 


V.  Form  of  the  output  data. 


The  results  of  the  computations  are  punched  on  a single  card  for  each 
source -sample  combination.  In  Appendix  E is  given  a listing  of  the  answer 
cards  resulting  from  the  computations  performed  with  the  data  shown  in 
Appendix  D.  As  in  Appendix  D,  the  columa  numbers  are  included  for  refer- 
ence as  the  first  two  lines  at  the  top  of  the  page.  They  are  not  punched 
out  in  the  normal  running  of  the  program. 

It  will  be  noted  in  Appendix  E that  the  data  are  located  in  the  first 
seven  fields  of  the  card.  The  data  are  punched  out  in  the  following  order : 

X,  first  field j I,  second  field;  Z,  third  field;  x,  fourth  field;  y,  fifth 
field;  z,  sixth  field;  and  identification,  seventh  field.  The  identifica- 
tion which  is  punched  in  the  se'venth  field  on  output  is  the  same  number  as 
punched  in  the  sec ond  field  of  the  control  card  on  input,  and  this  number 
therefore  correlates  the  results  with  the  proper  source-sample  combination. 

VI.  Suianary. 

A con^juter  program  has  been  described  which  will  allow  the  conversion 
of  spectral  directional  reflectance  or  spectral  transmittance  data  of  a 
sample  illuminated  by  any  light  source,  in  eluding  sources,  such  as  fluores- 
cent lanps,  which  contain  spectral  lines.  A detailed  description  of  how  the 
various  data  are  to  be  prepared  for  introduction  into  the  computer  has  been 
given  together  with  a description  of  the  output  data  from  the  computer.  The  i 
program  has  been  comple tely  checked  and  is  able  to  complete  each  source- 
sample  combination  in  an  average  time  of  approximately  1 minute  1$  seconds 
per  case.  Read- in  time  for  the  program  deck  is  approximately  $ minutes. 
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APPENDIX  A 


Listing  of  the  source  program  which 
was  assembled  by  the  FOR  TRANSIT  l(S) 
sygtem.  This  assembly  produced  the 
program  to  be  used  on  liie  IBM  6^0 
conputer . 
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SPECTRAL  LINE  COLORIMETRY'"'"' 

PROJECT  N0.60-4A9 
J.P. MENARD  650-CODE 

1 DIMENSIONA(40  ) *B(40)  ,C(40)  > 

IT (40 ) »AA( 5 ) »BB( 5 ) »CC ( 5 ) » 

_ 2TT ( 5 ) vWLLj  40)  »WL(  ) »WK( 5)  * 

3S ( 40  ) 

40  READ, A _ __  ^ 

45  RFAD»B  

50  PEAD,C 
55  READ,WLL 

60  READ  10,NN,M  

" 6 5 " GO  J 0""(  7 0 V 9 0 V 5 6 0 ) '»  N N 
_70  READ,S<r: 

75  READ  10, N 

80  READ,WL 
'“85  READ,WK 

90  READ  ; I 

Too'  j=r"  " - i vO 

105  DO  135  I=1,N  ' 

no  IF(WL(  I )"-WLCri  ) )560,115,'3d5 

115  AA(  I )=A(  J)  

120  BB(  I )=B(  J)  ' 

125  CC( I )=C( J)  ' r 
130  TT(1  )=T( J)"“  “ ■ ' 

135  CONTINUE 

140  EXX=0,0 

145  EYY=0.0 

150  EZZ=0.0 

155  EVV=0.0 

160  DO  T8  51= 1,40  " 

165  D = S(  I )*T(  I ) 

170  FXX  = FXX+A{I  )"*D 

175  E\  Y = EYY+B(_I  )*D 

180  FZZ  = FZ7+C(  n^D 
1R5  EVV  = EVV  + B(  I )*S(  I ) 
rod"  DO" 2 'l  51=1  , N 
195  D = WK ( I )*TT(  I ) 

200  EXX=EXX+AA (I ) *D 
205  EYY=EYY+BB( I ) *D 
'210  EZZ  = EZZ+CC(  I ) ^^D  ■ 

215  EVV  = EVV  + BB(  I )*'a/K  I I ) 

220  fx  = ex'x7fvv 

225  EY=EYY/EVV 

230  FZ  = FZZ/EVV  ' " 

235  D = EXX+EYY  + EZ7  

?40  EXX=FXX/D 

245  EYY  = EYY/D  _ _ 

250  EZZ=EZZ/D 
260  PUNCH  16,EX,EY,EZ, 
iexX,EYY,EZZ,M  ' 

300  GO  TO  60 

■305  ! F ( W L { r ) - W L L T2  )“  )“3“  10  ,"34  O' , "3  5 0 

310  J=1 

3 1 1 D = ( W LTTT-W  L 'LT  jTT7  ”■" 


- Ill  - 


1 (WLL  U+1  )~WLL7  J)  ) 


315 

320 

325 

AA(  I )=A(  J)+D*(A  ( J+1  )~A(  JM 

RR( I )=B( J+1 )-B( J)  ) 
CC(  I )=C(  J)+D*(C(  J+1  )“C:(  J)  ) 

330 

335 

TT( I )=T( J)+D*(T(J+1)“T( J)  ) 
GO  TO  135 

340 

345 

J=2 

GO  TO  115 

350 

355 

IFIWLI I )-WLL(40) )355t545*560 
IF( WU I )-WLL( 39 )) 360*525.535 

360 

365 

IF(WL( I )-WLL( J+1 ) >365*505.515 
D=(WL(  n-WLL(  J)  )/ 

l(WLUJ+l)-WLL(Jn 

370 

NN  = 0 , 

375 

•—1 

1 

< 

II 

5: 

< 

1 

380 

A0=A(J) 

385 

A1=A(J+1) 

390 

A2=A ( J+2 ) 

395 

GO  TO  600 

400 

AA( I )=AAA 

410 

AM1=B(J”1) 

415 

A0  = B ( J)^ 

420 

A1=R(J+1> 

425 

A2=B( J+2) 

430 

GO  TO  600 

435 

RB( I )=AAA 

440 

AM1=CIJ“1) 

445 

A0=C ( J) 

450 

M=C(  J+1) 

455 

A2=C(J+2) 

4673 

GO  TO  600 

465 

CC( I )=AAA 

470 

AM1=T(J-1) 

475 

A0=T ( J) 

480 

A1=T ( J+1 ) 

485 

A2=T ( J+2) 

490 

GO  TO  600 

495 

TT( I )=AAA 

500 

GO  TO  135, 

505 

J = J+1 

510 

GO  TO  115 

515 

J=sJ+l 

520 

GO  TO  360 

525 

J = 39 

530 

'GO  to  1-15 

535t 

J*39  • 

540 

GO  TO  311 

545 

J = 40 

550 

GO  TO  115 

560 

STOP 

600 

AAA=(D/2.0)*( A2-3.0 

:i*Al+3*0*A0-AMl  ) 

601 

AAA= ( AAA+( 2.0+AM1-5.0 

1*AO+4,0*A1-A2 ) /2.0 )*D 

602 

AAA= ( AAA+( Al-AMl )/2»0) 

- 15  - 


i*D+AO 

610  _NN  = NN  + 1 _ 

620  GO  TO  (466*435*46'5  *495  ) *NN 

FND__..,„_. 
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APPENDIX  B 

Listing  of  the  modified  source  program  which 
was  assembled  by  the  FORTRAN  II  system.  This 
assenbly  produced  a program  to  be  used  on  the 
IBM  70U  computer  in  simulation  of  the  IBM  6^0 
computer  and  served  as  a program  check. 


rk  n ri 


- 17  - 


SPECTRAL  LINE "COLORIMETRY 

PROJECT  N0.60-A49 
J.P. MENARD  650-CODE 
1 DIMENSIONA(40 ) *B(40) *C(40) * 
1T(40) *AA(5)*Be(5)fCC(5)  » “ 

2TT(5 ) ♦WLL(40) *WL ( 5 ) ♦ WK ( 5 ) » 

35(40)  “ “ 

10  FORMAT(2I10) 

" 12  F0RMAT(7F10.5 ) " ”"  ” 

16  FORMAT(6F10.5*I10|  ' 

40  READ  12»  ( A(  I ) ,"I  = 1»40) 

45  READ  12*(B(I) ,1=1,40) 

50"  READ — 12  ♦ red  rn  = 1T40  I*  

55  READ  12*(WLL( I )fl*l»40) 

60  READ  10»NN,M 

65  GO  TO  (70.90*560) fNNj 

70  READ  12.  (S(  I ) .1  = 1.40) 

75  READ  10. N S ' 

~ 80  READ  1 2 r ( WL  ( rr.  I = 1 .X  )"^  “ 

85  READ  12.(WK(I ).I=l.fn 
90  READ  12. fTrr).I=T,40) 

100  J=1 

105  DO  135  I=1.N 

no  IF(WL(  I )-WLL(  1 ) )^60. 115.305 

115  AA( I )=4(J) 

120  BB( I )=B( J) 

"125  CC(  I )=C(  J)  "h 

130  TT( I )=T( J)  ' ! 

135  CONTINUE  i 

140  EXX=0.0 

— 145-eyy=o;o — ^ 

150  EZ2=0,0 

““  155  EVV  = 0*0  ^ 

160  DO  1851=1.40 

~ - 165  D = S(  I )»T(  I ) 

170  EXX=EXX+A( I )*D 

17 5 "E  Y Y = E YY +B  ( I ) # D 

180  EZZ=EZZ+C( I )*D 

185  EW=EVV+B(  I )*S(T) 

190  DO  2151=1. N 

■'195  D=WK(  f)*TT(  I ) - 

200  EXX=EXX+AA( I )#D 

— 205EYY=EYY+BB( I )#D“ 

210  EZZ=EZZ+CC( I )*D  ? 

215  EVV=EVV+BB( I )*WK( I ) 

220  EX=EXX/EVV 

225  FY  = EYY/EW  ‘ 

230  EZ=EZZ/EVV 
■ 235  D=EXX+EYY+EZZ 
240  EXX=EXX/D 

245  EYY=EYY/D  

250  EZZ=EZZ/D 

“ 260  PUNCH  16.EX.EY.EZ. 

lEXX.EYY.EZZ.M 
-300  GO-TO  60 


305“I F ( WL ( n -WLLT2  y ) 3 1 0V34Q V 3 50 
310  J=1 

~311  Ds(WLn)-WLC(Jn/  ~ 

1 (WLL(  J+D-WLL  ( J)  ) 

315  AA  { I)  =A  ( J ) +D*  f A ( J+1 ) -A  ( JD 
320  BB(  n=B(  J»+D*(B(J+1  )-B(  J)  ) 

325  crii  r=C  ( JT+D*  (Ca^'Tr-^^ 

330  TT(  I )=T(  J)  ♦•D*(T(  J+l)-T(  J)  ) 

335  GO  TO  135 

340  J=2 

345  GO  TO  115 * 1 ^ 

350  IF(WL( I )-WLL(40) )355f545*560 
3 5 5~I  F ( W L ( I ) -W  Ctr(^9  )13  60  r52Tf 5 3 5~“ 
360  IF(WL( I )-WLL( J+1) ) 365*505 >515 

“365  D = ( WL  ( I ) -WLLIV) TA 

1(WLL(  J+D-WLU  J)  ) i 

'370  NN  = 0 ----‘-t 

375  AM1  = A( J-1 ) j ' 

380  “A0  = A ( J ) : 

385  A1=A{J+1) 

'390  A2  = A(  J+2  ) 

395  GO  TO  600, 

400  AA(I)sAAA 

410  AM1=B(J-1) 

-415  A0  = B(J) ; 

420  A1=B(J+1)  i 

425  A2  = B( J+2)  7 i“ 

430  GO  TO  600  ‘ ; 

435  BB(  I )=AAA : 

440  AM1=C(J-1) 

445-^0=C(^) 

450  A1=C(J+1) 

455  A2  = C(J+2) 

460  GO  TO  600 

465  CC(  n=AAA ^ 

470  AM1=T( J-1 ) 

475  "A0  = T (1) 

480  A1=T(J+1) 

485  A2  = T ( J+2  )“- 

490  GO  TO  600 

-495  TT(  I )=AAA 

500  GO  TO  135 

305^^  ^a+i 

510  GO  TO  115 

515  11= J+1 

520  GO  TO  360’ 

525  J=39  ^ 

530  GO  TO  115 

535~J=39  — 

540  GO  TO  311 

545  J=40 

550  GO  TO  115 

560  STOP ^ 

600  AAA= ( ( (D/2.0 )#(A2-3.0*Al+3.0* 

1 AO-AM 1 ) + ( 2T0* AM 1-510* AO+4 10*^ 
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lAl-A2)/2,0)»D+(Al“AMl)/2.0)* 

lO+AO 

'610  NN*NN+1 

620  GO  TO  (400*435»465*495) *NN 
' END(0*0*'1*1*0)  ---  


APPENDII  C 


A listing  of  the  data  constants 
required  by  -the  program. 
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APPENDIX  D 

A listing  of  typical  input  data  cards. 
(Note:  The  column  numbers  appearing 

as  the  first  two  lines  on  each  page 
are  not  a part  of  the  input  data.  They 
serve  in  this  appendix  only  to  index  the 
80  columns  available  on  each  card.) 
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APPENDIX  E 


A listing  of  answer  cards  resulting 
from  the  computations  using  the  data 
of  Appendix  D.  (Note:  The  column 
numbers  appearing  as  the  first  two 
lines  on  the  page  are  not  punched  out 
on  output.  They  serve  as  reference 
oily  in  this  appendix. 
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